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Abstract 

We have performed a detailed Monte Carlo study of a diffusionless (1 + l)-dimensional solid- 
on-solid model of particle deposition and evaporation that not only tunes the roughness of an 
equilibrium surface but also demonstrates the need for more than two exponents to characterize 
it. The tunable parameter, denoted by ^, in this model is the dimensionless surface tension per 
unit length. For /i < 0, the surface becomes increasingly spikier and its average width grows 
linearly with time; for = 0, its width grows as \/t. On the other hand, for positive /z, the surface 
width shows the standard scaling behavior, (crm(i)) ^ f{t/M^/^) where M is the substrate 
size and f{x) const {x^) for x large (small). The roughness exponent, a = 1/2 for fi < 2, and 
= 3/5,4/5 & ~ 1 for /X = 5, 6 & 7 respectively; the growth exponent, f3 = 1/4 for /i < 2 and 
= 1/2 for fi >~ 3.5 respectively. These exponents are different from those of the height-difference 
correlation function, q' = 1/2, /3' = 1/4 and z' = 2, for higher values of fi suggesting thereby that 
the surface could be self-constraining. 
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I. INTRODUCTION 



Stochastic evolution of rough surfaces into stationary self-affine structures has been a subject 
of intense study . The initial growth and the asymptotic stationary value of the width, 

{(^mii)), of a surface are described by the dynamical scaling form, ~ M"/(t/M"/^) 

where M is the substrate size and /(x) const (x^) for x large (small). The scaling 
exponents /5 and a lead to a classification of the surfaces into distinct universality classes 
such as, for example, the Kardar-Parisi- Zhang (KPZ) class ^ and the Edwards- Wilkinson 
(EW) class 

The simplest physical process giving rise to the growth of rough surfaces is the Random 
Deposition (RD) of finite sized solid particles on a fiat substrate; at every instant of time, the 
height of the surface at a randomly chosen site is incremented by the size of the particle - i.e., 
a Solid-on-Solid (SOS) event occurs at a random site. Since the SOS events are statistically 
independent of each other, the surface growth is equivalent to a one dimensional random 
walk for the height variable; consequently, the surface- width, {amit)), is independent of the 
substrate (or system) size and increases monotonically with time, {am{t)) ~ t^^^. 
On the other hand, if the particle sticks to the surface as soon as it makes the first contact, 
then the direction of growth is not always perpendicular to the substrate as it is in the case 
of RD process. Possibility of such lateral growth leads to the formation of voids in the bulk 
and overhangs in the growing surface. The width of the surface, however, saturates to a 
(substrate) size-dependent value because of a build up of spatial correlation between the 
surface heights. Surface grown by such a process, known as the Ballistic Deposition (BD), 
belongs to the KPZ class. 

An SOS event followed by the diffusion of the deposited particle on the surface smoothens 
the surface. The morphology of such smoothened surfaces depends on the nature of diffusive 
transport taking^jlace - namely, whether it is due to local height differences jsl or due to 
local curvature j6|. Most of these models use simple phenomenological or adhoc rules for 
the diffusion of deposited particles on the surface. The asymptotically saturated surface is 
referred to as the equilibrium surface, not in the thermodynamic sense but in the sense that 
its width has attained a stationary value. 

The growth of such a surface can also be described by the height-difference correlation 
function (HDCF), GM{k,t), that is defined in terms of the average differences in the heights 
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of sites separated by a distance 0,0] - namely, GM{k, t) = {[h{i + k, t) — h{i, t)Y) where the 
averaging is over all reference sites i and also over independent realizations of the surface. It 
has been shown |7|,l8| that GM{k,t) ~ k'^''' fcik/^Git)), where a' is known as the wandering 
exponent; initially {t <^ M^'), the correlation length, ^aii)^ increases with time as ^oii) ~ 
t^/^' but saturates to a value proportional to M for t ^ M^'. Normally, for SOS models with 
local growth rules, a' = a and z' = z. However, for SOS models with global constraints on 
surface heights such as, for example, the Suppressed Restricted Curvature model [9|, a' 7^ a 
and z' 7^ z, eventhough their ratios may be equal {a'/z' = (3 = a/ z). An interesting question 
to ask is whether there is a growth process that necessitates the need for more than two 
exponents to characterize an equilibrium surface. 

In this paper, we present a detailed Monte Carlo study of a diffusionless (1 + l)-dimensional 
model of particle deposition and evaporation that not only generates an equilibrium surface 
with tunable roughness but also demonstrates the need for more than two exponents to 
characterize it. This model is based on the observation that adding or removing a particle 
is equivalent to changing the area of the surface by one basic unit. It has a single tunable 
parameter, denoted by /i, the magnitude of which measures the strength, in units of thermal 
energy, of this change. 

We show that the morphology of the one dimensional surface changes drastically from spiky 
to rough as /i, the surface tension parameter, changes sign. More interestingly, we show 
that different ranges of /i values correspond to different universality classes - (i) EW class 
for yU in the range < /i < 2; (ii) growth exponent (5 = 1/2 for /i > 3.5 while the roughness 
exponent a increases monotonically with /i - namely, a ~ 1/2, 3/5, 4/5 and 1 for yU = 3.5, 5, 6 
and 7 respectively. The Monte Carlo data presented here suggest that the parameter n can 
tune the exponents a and (3. On the other hand, the exponents a', z' and /?'(= a'/z') of the 
HDCF, GM{k,t), retain their EW values independent of fi. They are, thus, different from 
those of the surface width eventhough no arbitrary global constraint is imposed on surface 
heights. This is in contrast to the Suppressed Restricted Curvature model 9|] that explicitly 
invokes a global constraint on the extremal height of the surface. The paper is organized as 
follows. 

In section II, we introduce and discuss the model in detail; in section III, we present a detailed 
Monte Carlo study of the one dimensional model and finally, in section IV, summarize the 
results. 
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II. A (1 + l)-DIMENSIONAL SOS DEPOSITION-EVAPORATION MODEL. 



Let A4{t) denote a surface configuration at time t growing on a one dimensional lattice 
substrate consisting of M sites. The particle being deposited (or removed) is assumed to 
be a unit square so that the surface Ai{t) is actually a collection of integer- valued variables 
(heights), {hj{t) \ j = 1,2, ■ ■ ■ , M}. The contour length of A4{t) can then be written as, 

M 

LM{t) = M + J2\hi+iit)-hi{t)\ (2.1) 

i=l 

with periodic boundary codition, hM+iif) = hi{t), invoked. 

During the growth process, if a particle is deposited or removed at a randomly chosen site 
j, then the change in the contour length because of a change in hj{t) by ±1 is given by 

AL^(t) = L'j^it) - LmH) (2.2) 
= I /i,-+i(t) -/i;.(t) I + I h'^it) - hj_i[t) I 

- I h,+i{t) - h,{t) I - I hj{t) + hj_,{t) I (2.3) 

where h'j{t) = hj{t — 1) ± 1, say with probability b or 1 — b respectively; periodic boundary 
condition may be imposed for the edge sites 1 and M. It is clear from this definition that 
AL{t) e {2,0,-2}. 

Acceptance of the new value, h'j, of height at site j could be on the basis of the standard 
Metropolis criterion - namely, that it is accepted with probability p[hj] = min(l, e"'^^'^^*-'), 
where /i is a tunable parameter in units of ksT. The surface evolves stochastically as the 
sites are scanned and a choice for random deposition or removal of a particle is exercised. 
The excess contour length, L^it) — M, defined in Eq.(l) may be recognized as the energy 
of a solid-on-solid (SOS) model [lO|, 

E = fxJ2\Ht)- hj{t) 1= ^l{LM{t) - M) (2.4) 

<ij> 

with which the surface configuration M. is weighted. It is clear that /i may be interpreted 
as a microscopic surface energy ("surface tension"), and the summation is over nearest 
neighbor pairs. In fact, it was used 11, [l2| for studying the influence of an external field on 
the equilibrium fluctuations of an overhangs-free interface in the Ising model. This (1 + 1)- 
dimensional model may also be considered as that of a directed polymer with fluctuating 
number of monomers. 
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FIG. 1: Surface profiles after 20000 Monte Carlo sweeps for fi 
top to bottom). Substrate size M = 128 
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Typical surface profiles obtained after 20000 Monte Carlo sweeps are shown in Fig.([T]) for /i 
in the range —0.1 to 2.0. They are smoother for positive than for negative fi. In fact, the 
surface becomes spikier as fi becomes more negative. In all these cases, the mean height 
of the surface, (hm), does not change because deposition and evaporation are equally likely 
events. 

For a given realization of the surface, the width (or equivalently, the roughness) is given by 
the mean square fluctuation in the surface heights, 

M 



(2.5) 



where h(t) denotes the center of mass, or equivalently, the mean height of the surface: 



1 



(2.6) 



Averages over a number of realizations may be denoted by (cr|j(t)) and (/iA/(t)) respectively. 
Since the model is defined in terms of the contour length, L, temporal behavior of the 
fluctuations in L should be of interest. However, it is not difficult to see that they are 
directly related to those of surface heights: 

Let us denote the deviation of a column height from the mean height (center of mass) of the 
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surface by Hj{t) = hj{t) — h{t). Then we have, from Eq. fl2.ip 

1 



^{mt)f) = mt)-Mf) (2.7) 

M 



= Y,m-.i{t) - K-,{t)f) + m (2.8) 
1=1 

= 2M{aUm-C,{t)]+X(t) (2.9) 

where 

/ M /M-i+l \ \ 

^(^) = 2 I H,{t) - H,^,{t) I 5^ I i7,+i(t) - H,{t) I j ^ (2.10) 

consistes of terms proportional to the correlation functions, Ci(t), C2(t), ■ ■ ■ , and hence can 
be written as 

X{t) = 2M{al{t)) /[Ci(t),C2(t),C3(t)---] (2.11) 
The correlation functions, Ci(t), C2(t), ■ ■ ■ , are given by the standard definition. 

In the case of a pure random deposition-evaporation process, the individual columns are 
statistically uncorrelated with each other and so, we have the simple result, 

l{{5L{t)y)^{al,{t))^t (2.13) 

In all the other cases (/i 7^ 0), the time-development of the correlation functions cannot be 
ignored. 

III. RESULTS AND DISCUSSIONS 

Starting from the initial straight line configuration, we have monitored the build-up of the 
surface roughness by estimating the time-dependence of the root mean squared fluctuation, 
(o'm(^)) = (^m(^))^^^) averaged over many independent runs employing periodic boundary 
conditions. 

We have presented in Fig. ([2]) the Monte Carlo estimates of {(Tm{t)) for both positive and 
negative values of fi. It is clear that {(Jmit)) is proportional to t for /i < 0, whereas it is 
proportional to t^^^ for > 0. This suggests that the surface may belong to the EW class 
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for positive fi, in which case {am{t)) will have to saturate asymptotically to a value that 
scales with the system size M. 

Following previous works [isl, the finite size scaling of (cTm(t)) for > is given by 

f : t < 

{am{t)) = f{tM-') ; /? = a/z, (3.1) 

where f{x) ~ z = const for a: 3> 1 and /(x) ~ x"/^ for x -C 1. 



A. Dynamics of roughening (0 < /x < 2) 



The scaling analysis of our simulation results, presented in the form of data collapse for 
fi = 0.5 and 2.0 in Fig.Q, clearly indicate that the exponents have the values expected of 
the Edwards- Wilkinson (EW) class 0, [s]. 
z = 2. 



13| for < /i < 2 - namely, (3 = 1/4, a = 1/2, and 
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FIG. 2: (a) Time-development of root mean-square fluctuation, am{t), for /i = 
-1.0, -0.75, -0.5, -0.25 and 0.0 (top to bottom). System size M = 128. Slope of the line, /3 = 1/2 
for fi = 0, and is asymptotically equal to unity for the other negative values of /i. (b) am{t) for 
positive values of = 0.0,0.1,0.03125,0.0625,0.125,0.25 and 0.5 (top to bottom). System size 
M = 512. Clearly /? = 1/4 for /i > 0. 
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FIG. 3: Data Collapse for root mean squared height-fluctuation, (a) fx = 0.5 and system sizes 
M = 32 to 128 in steps of 16. (b) fj, = 2 and system sizes M = 64 to 256 in steps of 32. The 
exponents are a = 1/2 and /? = 1/4 in either case. 



The /i-dependence of the saturated value, {(Jm)s, is reflected in the proportionality factor, 

{am)s = B{fi)M'' (3.2) 

where a = 1/2 for < < 2. Numerical evidence for this is presented in Fig.(jl^) in which 
size dependence of {o'm)s is shown for various values of /i. Moreover, data presented in 
Fig. (jib) for system size M = 128 suggest that B{fi) is an exponentially decaying function, 

5(/i) = 5o e-^/^^ (3.3) 

where the constants are Bq = e^-^^^'^^ ''^ and jjB = 1-22 ± 0.03. 

Since fluctuation of contour length is linearly related to that of surface heights (Eq. fl2.9l) ). 
it will be of interest to know the dependence of the average contour length, (L), on system 
size and fi. Data presented in Fig.(l5H) for /i = 0.1 and 2 suggest that (L) is proportional to 
the system size: 

(L) = A{^i)M. (3.4) 

The prefactor, A{fi), is an exponentially decaying function of /i, as is evident from Fig.(l5b) 
for system size M = 128: 

A(/.) = e-^/'^^ (3.5) 
8 



where the constants are Aq = 6.49 ± 0.13 and fiA = 3.2 ± 0.2. 




(a) (b) 

FIG. 4: (a) Average saturation values, {(Tm)s as functions of sytstem size M for ^ = 0.25, 0.75 
and 2.0. Slope of the fitted straight lines is 1/2 within errors, (b) Dependence of the equilibrium 
values {(Tm)s on ^ for a system size M = 128. 



B. Dynamics of roughening (/i > 2) 

Interestingly, we find that the exponents, a and /3, do not retain their EW values (1/2 and 
1/4 respectively) as yU is increased beyond the value 2. This may not be surprising for the 
following reason. 

Since the probability of growth, and hence of increasing the contour length, at any given 
site is proportional to e~^^ (Section II), larger the value of yU lesser will be the probability 
for growth. In the extreme limit fi ^ oo, the surface hardly evolves from its intial flat 
morphology and at best may acquire a few steps in the long time limit t — > oo. In other 
words, it will remain mostly fiat, save for a few steps here and there, and so will not resemble 
an EW surface. It is therefore reasonable to expect that a will change from its EW value 
to that corresponding to a surface with fiat morphology as we increase /i. 
In Fig. ([6]), we have presented data collapse for /i = 3.5, 5, 6 and 7. We see that the collapse 
obtained is reasonably good for system sizes ranging from M = 128 to 512, though the 
statistics becomes poorer especially for small system sizes at higher values of fj,; because of 
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FIG. 5: (a) Size-dependence of the saturation values of the average contour length {L)s for /_f = 0.1 
and fi = 2.0. In either case (L) oc M; the ;U-dependence is only in the proportionality constant, (b) 
Semi- log plot of the average contour length {L)s for M = 128 fixed as function of surface tension 
/i. 

small e~^'^, larger systems and longer time spans are needed for acceptable statistics on the 
build-up of fluctuations. That is why, collapse is shown Fig. ([611) for only three values of M 
( = 512,448 and 384) at /i = 7. 

The growth exponent /9 = 1/2 for all these values of /i, whereas the roughness exponent a 
increases monotonically with /i (a ~ 1/2, 0.6, 0.8 and 1 for fi = 3.5, 5, 6 and 7 respectively). 
Larger the value of /i, lesser will be the probability for a change in the surface height that 
leads to an increase in the contour length [i.e., surface area in one dimension); consequently, 
we expect the growth of height fluctuations to be slow. Yet, /3 has increased from its EW 
value (= 1/4 for /i < 2) to that corresponding to Random Deposition (= 1/2 for fi > 3.5). 
In this context, it will be of interest to look at the height-height correlation data as well. 
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FIG. 6: Data collapse for various values of M ( = 512, 448, 384, 320, 256, 192 and 128) (a) fi = 3.5, 
and the exponent values are (3 = 1/2 and a = 1/2. (b) /i = 5, and the exponent values are (3=1/2 
and a ~ 0.6. (c) /x = 6, and the exponent values are (3 = 1/2 and a ~ 0.8. (d) fi = 7; data 
presented for M = 512, 448 and 384 = 5. The exponent values are (3 = 1/2 and a ~ 1; 



C. Height-height correlations. 

The height-difference correlation function (HDCF) 0, Q] of a surface configuration at time 
t - defined by the set of heights, {hi{t) \ i = 1,2, ■ ■ ■ , M} - is given by 

GM{k,t) = -J2{ih,Mt) - h,it)f) (3.6) 
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It shows a scaling behavior, 



Guik, t) ~ k^'-'fc ( -r^TT ) ; A; = 1, 2, ■ ■ ■ M/2 (3.7) 



Xcit) 

with respect to the lag variable k. The exponent a' is known as the "wandering" exponent. 
Normally, after growing as t^^'^' at early times, the correlation length ^g(^) saturates to a 
value proportional to the system size: 

f t^/'' for t^/^' < M 
CG{M,t)^{ (3.8) 
y M for t^/' > M. 

For SOS models with local growth rules, a' = a and 2;' = z. 

It is clear from Eqs. fl3.6[ 12. 5p that GM{k,t) and {cj^it)) are related to each other by the 
following exact identity: 

M 

M 



GM{k,t) = l^Y.{[h,^,(t) - h,{t)f) 
1 ^ 



M 
i=i 

= 2((al,(t))-CM(A;,t)); A; = 1,2,---,M (3.9) 

where hi{t) = hi(t) — hi(t) and Cuik^t) is the height-fluctuation correlation function 
(HFCF), given by the standard definition, 

1 ^ 

CM{k,t) = —Y^{h,+,{t%{t))- CM(0,t) = {alM (3.10) 
i=o 

From Eq. (l3.6p . it follows that 

GM{Q.t) = = GM(M,t) (3.11) 

Gm{M/2 -k,t) = Gm{M/2 + k,t); A; = 1, 2, ■ ■ ■ , M/2 - 1 (3.12) 

which implies that GM{k,t) is also symmetric about M/2, namely, 

Cm{M/2 -k,t)= Gm{M/2 + k, t) (3.13) 

Moreover, summing GMik,t) over /c, we get 

1 

-EGM(A:,t) = 2(a^,(t)) (3.14) 



M 
fc=i 
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which, together with the identity Eq. fl3.9p . imphes that 

1 

-J2CM{k,t) = (3.15) 

k=l 

Monte Carlo averages of CM{k,t), obtained after a time t = e^^'^ and normahzed with 
respect to Cjvf(0,t) = {a\^{t)), are presented in Fig. ([7]) for a system of size M = 512 and for 
fi = 1/2 & 5. It is quite evident that CM{k,t), is an exponentially decaying function over a 
reasonable range of k: 

CM{k,t) = Ace-"/^^^''''^ (3.16) 

where the proportionality constant Ac will depend on /i while the correlation length, 
^c{M, t), will depend on M and t. Monte Carlo estimates of the correlation length, ^c{M, t), 
are presented, and compared with those of C,g{M, t) in Fig.([8K,b) for = 1/2 & 5. It is clear 
that 

fortVs^M 

^c{M,t)<x^G{M,t)^ { (3.17) 

I M for ^ M. 

In view of Eq.f l3.8p . this implies that z' = 2. 

The exponential decay of the HFCF, CM^k, t), from a positive value {cr\i{t)) at = 0, taken 
with the equality Eq. fl3.15l) . implies that CM{k,t) has to be negative over a certain range 
k G [M/2 — ko, M/2 + ko] for some k^. In turn, this implies (by Eq. (13.91) ) that Gjuik, t) has 
to be an increasing function of k. In fact, GMik,t) has been shown to have the following 
scaling behavior [7, 



For t > M^ 



, for A; < M 

GM{k,t)r^{ , (3.18) 

for k ^ M/2. 



and for t < 



, P"' forA;<ti/^' 
G^M(A;,t)~<^ (3.19) 
forA;>ti/^. 

Monte Carlo estimates of GM{k,t) at various times for a system of size M = 512 and for 
yU = 1/2 & 5 are presented in Fig.([n^) in the form of a collapse diagram suggested by the 
above scaling forms. The corresponding saturation values are shown in Fig.(PD). Within 
numerical errors, we have a' = 1/2, j3' = 1/4 and z' = 2. It must be noted that these values 
are independent of /i, in striking contrast to those of the surface width. For example, for 
^ = 5,a^ 0.6, (3 = 1/2 and z ^ 1.2. 
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In order to show that this is not a numerical artifact, we consider Eg. (13.90 in the initial 
time regime t -C M^. For small lag distances (say k ^ M/4), the HFCF CM{k,t) is an 
exponentially decaying function, Eg. (13. 161) . whatever be the value of fi. At a given time 
t, the correlation length ^c{M,t) oc t^^^' in the time regime of interest, and (o"|j(t)) is a 
constant that does not depend on k. We therefore have, for k <^ t^^^' , 



where Kc is a proportionality constant. Comparing it with the expected scaling behavior, 
Eg. (13. 191) . we have the result a' = 1/2 independent of /i. 



lapse of time t. Correspondingly, we see that Cuik, t) is more or less saturated to a negative 
value, say -C^^]{t). The saturated values, estimated as averages of the last M/8 data points 
at various times, for /i = 1/2 & 5 and M = 512 are presented in Fig. (fTUk .b). Straight line 
fits to the data lead to the estimate j3' ^ 1/4. The corresponding estimates of (o"|f(t)) are 
also presented in these figures. 

Since in this time regime, (o"|/(t)) ~ t"^^, we can write 



Which exponent will characterize the time-dependence of G^{t) is decided by the propor- 
tionality constants K„ and Kq- 

In the case /i = 1/2, straight line fit to the values of {<j\j{t)) in Fig. lfTOk ) give an estimate 
/5 ^ /5' ~ 1/4. Hence, G^^]{t) oc t^^^ - independently supported by the Monte Carlo data as 
well. 

On the other hand, for /i = 5, /3~l/2 and /?' ~ 1/4 are the estimates obtained by fitting the 
data presented in Fig. ffTOb ). The proportionality constants have the values, K^j ^ 2.5 x 10~^ 
and Kc ~ 0.005, respectively. From Eg. (13.211) . it is clear that G^^]{t) ~ t"^^ only for times 
t ^ {Kc/Ka)"^ ^ 400. However, data presented in Fig.(PD) indicate that {cr\[{t)) saturates 
for t ^ 400. Hence, G^lj{t) oc t^^^ - independently supported by the Monte Carlo data as 
well. 




(3.20) 



For lag distances k ^ t^/^'. 



we observe that GM^k, t) saturates to a value, say G^/{t), after a 



lG'^^{t)=K^t'^ + Kct'^' 



(3.21) 
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FIG. 7: Semi-log plot of the Height-Fluctuation correlation function, Cjvf(A;,t), for M = 512 and 
/i = 1/2 & 5. The data were obtained after a time t = e^^'^ Monte Carlo sweeps, and averaged 
over 64 independent runs. It is clear that CM{k,t) is an exponentially decaying function over a 
reasonable range of lag distances {k G [0,~ 48]). 

D. Discussion. 

Metropolis dynamics of surface growth based on the SOS Hamiltonian, Eq J2.4l has aheady 
been studied by Siegert and Phschke In their model, particle deposition is a random 
event that is extraneous to the energetics of the evolving surface, whereas diffusion of ran- 
domly deposited particles is what relaxes the surface towards an equilibrium energy state. 
Their model is motivated by the observation that relaxation of a surface, being grown by the 
molecular beam epitaxy process, is through the diffusion of particles on the surface. They 
have shown that the surface belongs to the EW class. 

On the other hand, in the model we have presented here, the surface relaxes towards its equi- 
librium energy state not by the process of diffusion but by the process of particle deposition 
or evaporation. Since a change in energy is associated with an event of particle deposition 
or evaporation, we may say that this model treats a growing surface as a canonical object 
evolving towards its equilibrium state. Again, since equilibrium energy is the same as the 
equilibrium contour length (surface area, in higher dimensions) which, in turn, is equivalent 
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(a) (b) 

FIG. 8: (a) Correlation lengths, £,G,c{M,t), for a system of size M = 512 as a function of time. 
Filled circles correspond to ^G{M,t), while filled triangles correspond to S,c{M,t); upper ones are 
for fi = 5 and the lower ones for = 1/2. Within errors, S^G,ciM,t) ~ t^/^. (b) Correlation lengths, 
^G,c{M,t), as a function of M after a time t = e^^'^ Monte Carlo sweeps. Filled circles correspond 
to ^G{M,t), while filled triangles correspond to ^c'(-^;0; upper ones are for = 5 and the lower 
ones for fi = 1/2. Within errors, S,G,ciM,t) ~ M. 



to the equilibrium number of particles defining the surface, we may also say that this model 
treats a growing surface ClS Sb grand canonical object. 

The motivation for this model is to provide a microscopic view of an equilibrium elastic 
membrane consisting of a number of folds of varying sizes for a given parameter fi (in 
thermal units). Therefore, diffusion cannot be the mechanism for surface relaxation in the 
SOS picture adopted here. The Monte Carlo study presented here provides a strong evidence 
that the parameter /i can tune the exponents, a and /3, characterizing the scaling behavior 
of the surface-width. 

Normally, in small scale simulations (M not very large), better estimates of these exponents 
are obtained from the scaling behavior of the height-difference correlation function (HDCF), 
GM{k, t), rather than from that of the surface- width. Jeong and Kim 9j have recently shown 
that this is true only for models such as the Restricted Solid-on-Solid (RSOS) model 0, Q] 
that have no global constraints on the column heights. Using a restricted curvature model, 
they have demonstrated that the exponents of GM{k,t) are not the same as those of the 
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Ln(k / 1"') Ln(t) 

(a) (b) 

FIG. 9: (a) Scaling collapse of the HDCF, GM{k,t), estimated at times ln{t) = 5,6,7, 8 & 9 for a 
system of size M = 512 and for /i = l/2&:5. It is evident that a' = 1/2 and z' = 2 independent 
of fi. (b) Saturation values of GM{k-,t) at these times for fi = 1/2 & 5. Straight line fits give 
/3' = 1/4, within errors, independent of //. 



surface width when the extremal height of the surface is suppressed. 

Eventhough no such global costraints are there in our model, the exponents of (^^/(/c,^) 
are different from those of the surface-width for yU > 2; the surface is self-constraining for 
higher values of /i. The subtle point this study has brought out is the decisive influence of 
proportionality constants on the scaling behavior of GMik,t). 

For example, the time-dependence of the saturation value of GM{k,t) (typically, for k > 
M/4), as given by Eg. (13.211) . depends crucially on the proportionality constants K„ and 
Kc. If /5 > then for times t > {Kc/K^y/^^f^-'^'\ G£^(t) oc t^f^ provided the surface 
width has not yet saturated; else, G^^]{t) oc t"^^' . In our model, f3' = 1/4 independent of n, 
whereas /5 = 1/2 for fi >~ 3.5. 

Since the height-fluctuation correlation function (HFCF) CM{k,t) is an exponentially de- 
caying function of the lag-distance k (Eg. (13.161) ). whatever be the value of /i, the eguality 
Eg. (13.91) for GM{k,t) can be rewritten as 

GMik,t)^{alit))-{l-^^^^y, k<^UM,t) (3.22) 

In the growth regime (t <C M^), at a particular time t, the fc-dependence of Guik^t) is 
given by GM{k,t) oc k; this implies, by definition Eg. (l3.19p . that a' = 1/2, whatever be the 
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FIG. 10: Saturation values of the correlation functions, GM{k ^ t^^^,t), CM{k ^ t^^^,t), at 
various times t{<^ M^) and the corresponding mean squared widths, (o'|^(t)). (a) fj, = 1/2: The 
proportionality constants are Kq ~ 0.415, Kc ~ 0.001 and K^j ~ 1.33. The exponents are 
P f3' 1/4. (b) fj, = 5: The proportionality constants are Kq ~ 2.11, Kc ~ 0.005 and 
~ 0.00025. The exponents are /3 « 1/2, /3' w 1/4. 



value of /i. In contrast, a increases with /i for /i >~ 4. 
IV. SUMMARY 

We have presented an SOS model of particle deposition and evaporation for describing the 
roughness of a (1 + l)-dimensional equilibrium elastic membrane with folds of varying sizes. 
The stochastic evolution of an intial flat surface towards its equilibrium configuration is 
governed by the standard Metropolis criterion based on the energy (or equivalently, the 
contour length), given by Eq. (12.41) . associated with this surface. The only tunable parameter 
in this model is fi, the dimensionless surface tension parameter. 

For negative values of fi, the surface becomes increasingly spiky and its width (root mean 
squared fluctuation of its heights), a oc t. But for /i = 0, a oc \/t. In either case, a does not 
saturate to a stationary value. 

For /i > 0, on the other hand, we observe an asymptotic saturation of a indicating that the 
surface evolves into an equilibrium state. In fact, our Monte Carlo data for < /i < 2 show 
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that the surface belongs to the Edwards- Wilkinson class characterized by the exponents, 
(3=1/4: and a = 1/2. However, for /i = 3.5, the growth exponent (3 is different {(3 = 1/2) 
while the roughness exponent, a, has the same Edwards- Wilkinson value. For /i >~ 4, we 
observe that a increases with /i, while (3 (= 1/2) remains constant. Our Monte Carlo study, 
thus, provides a strong evidence that the roughness exponent, a, of the (1 + l)-dimensional 
surface can be tuned by the surface tension parameter, /i. 

We have also shown that the exponents of the height-difference correlation function differ 
from those of the surface width for higher values of fi suggesting thereby that the surface 
could be self-constraining. We have presented a brief heuristic argument to support this 
observation. There is a definite need for a detailed analytical, and large scale numerical, 
study of this simple but rich model 
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